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ABSTRACT 



Theories of rough surface scatter and gas bubble behavior are used 
with the Pierson-Moskowitz wind-wave spectrun-i and an empirically- guided 
formulation of bubble concentrations at sea to calculate the true velocity 
gradient and losses "at the surface". These values are then entered into 
Bucker's wave theory solution for sound propagation and leakage in a 
surface duct. Curves of propagation loss are calculated for comparison 
with ocean test data obtained with the SQS-26 sonar. The predictions axe 
shavn to be significantly better than those based on the eirpirical equations 
of project AMDS. 



I. INTRODUCTION 



The purpose of this paper is to show that physical descriptions of 
bubble behavior arid sound scatter at the sea surface, extended by laboratory 
e:q>eriraents and in-situ measurements of near-surface bubbles have, at last, 
made it possible to give reasonably accurate predictions of surface duct 
propagation from a knowledge of wind speed over the surface and duct velocity 
profile . 

We will start with a review of the current oceanographic description 
of the sea (Section II), and surface scattering theory (Section III), supple- 
menting the basic formalism by the results of recent laboratory studies. The 
theory of sound dispersion, scatter and absorption by bubbles, simple bubble 
kinetics, the extrapolated results of in-situ near-surface bubble measurements 
are developed in Section IV. 

(U) Finally (Section V), and with the very generous cooperation of Dr. 

Homer R. Bucker of the Naval Research Laboratory, who used his computerized 
wave guide solution^ as a test for our surface loss and refraction predictions, 
we look again at the data of two SQS-26 propagation experiments performed by 
the Key West Test and Evaluation Detachment. Bucker 's wave theory solution • 
is modified by using the observed wind speed at the time of the experiment 
to correct the velocity gradient for bubble presence and to determine the 
bubble and surface losses which now replace an empirical one db surface loss 
factor that he used to obtain excellent data fits. 
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II. THE SEA 



We assume that the sea surface from which the sound is scattered is 
homogeneous, and that the temporal variation of sea heights, J , at a point 
is Gaussian distributed with mean value zero, 

<C> = 0 

2 

• and with variance , a , 



<C^> = 



The surface height correlation function is 



C = -2 



(C (0.0,0) C (u,v,T) > 



( 1 ) 



where u, v, “T are the spatial and temporal lags. 

The Fourier transform of 2 yields the three-dimensional spectrum of the sea 



X (x 



X 



X , f2) = 

y 



(2:r) 
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Z (u,v,t) exp 



Ci (x u 

X 



x^v - n t) ] dudvdT (2) 



where ~ components of surface wave propagation constant 

X y 

= angular frequency of surface wave 

We shall assume the gravity wave relation between angular frequency, propagation 
constant emd acceleration of gravity, g: 

= xg . 
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m 



■liWi 



This permits the directional spectrum of the sea to be specified either 
in the two-dimensional rectangular coordinate form 

X = >1^ (X ,a) 6 (vp-- n) (3a) 

or the polar form, 

•X = >1^ 6(V?X - Q) 

where the- Dirac delta function,^ (v'gx - fi) , simply selects the wave that 
fulfills the gravity wave relation betv;een frequency and propagation constant. 

The more commonly measured frequency spectrum is now obtained by 

using the energy density coordinate transformation 

/ ^]^dKda = dn 

0 



This permits the frequency spectrum to be expressed as the integral of the 
polar form of the X spectruin over all asinuth angles. 




Ah) 



The semi-empirical form of the frequency spectrum for the fully- 

Z 

developed sea that is most widely accepted today (Pierson-Moskowitz) 

is 




. ( 5 ) 



% 



U 



where 



a = 8. 1 X 10 







W 



19.5 



g 



0.74 



g/W 



19.5 



Wind speed at 19. 5M above sea surface 
acceleration of gravity in con-sistent units 



The three statistical parameters of the sea surface that are most significant 
in specifying the sound scatter are the mean square suface height, the mean 
square surface slope, t and the surface height correlation function, C. 



For a locally- generated sea we obtain the mean square height by 




For the mean square slope, it is convenient to use the relation 

3 

determined from optical measurements of the sea 



= 5.12 X 10”^W^^ + 0.003 + 0.004 



(7) 



where W^^^is measured in cm/sec at a height of 1^1 feet over the surface. 
Unfortunately the measured speed is sensitive to the height of the anemometerj 

h 

however corrections can be made. 
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Finally, for simplicity and because of the gener^ lack of two- 
dimensional sea correlation data, we will assume that the spatial correlation function 
is of Gaussian form and isotropic: 

C(jD = (8) 

where V- 

0 JT"7 ^ 

X = 'u + V 

and L = surface correlation length. 

The assumption of an isotropic Gaussian function makes it possible to calculate 
the correlation length from a knowledge of Titis slope, y . > rms height, or , 

by the relation that is valid for seas of lesser roughness : 

E = - 

L 



. ( 9 ) 
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ni. SOUIJD SCATTERING THEORY AliD LABORATORY STUDIES 






The theoretical solution to the sea surface sound scattering problem 
usually starts with the Helmholtz Integral; 



= i r 

L ^2 



— 1 

(D 

rr 

to 


_ 




L ''2 


Bn 


Sn V 1^2 /_ 



dS 



( 10 ) 



where p. 



k = 
X = 

H = 

n = 
See Figure 1, 



the scattered acoustic pressure at an interior field point, Q, 
of an enclosed volume, 

integrating surface selected so that acoustic field is zero 
everywhere except on the ensonified water surface. 

2^r/ X = acoustic propagation constant 
acoustic wave length in water 

distance from surface scattering region to field position, 
normal to surface scattering point. 



After applying the boundary conditions for acoustic pressure and 
toustic particle velocity at the sea surface, and assuming that there is no 
ladowing or secondary scattering even for near-grazing incidence, we 
Jtain the integral for the scattered pressure (now called A ?2 = Pq) from a given 

trface of scattering area A A in terms of the geometry: 



= Apj 



X. 

AA 



r --^dS 

plane ““ '^AA 



(11) 



nere Ap, 



j = scattered pressure for a plane surface AA over which 

‘plane there is negligible eorvature of the inciaent wave front. 

= Ikp^ AA cos e^^^'^2 ^i^V2:rr2 
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AA * -jmsonified surface area over which is constant and the 

sagitta of the wave front is less than A /8 • 






F 



o /incident sound pressure at sea surface. 

= incident sound angular frequency 

1 + cos 6^ cos 02 - sin 0^ sin 02 cos 0 
** cos 0^ (cos 0^ + cos Q^) 



e. 



0 , 



B. 






~2 






= angle of incidence (measured with normal) 

= angle of scatter (measured with normal) 

= azimuthal angle of scatter plane with respect to 
incident plane. 

= Ic (ij sin 0^ - ^2 cos 0^) 

= k (X^ sin 02 cos 83 + 12 ®2 ^ ®2^ 

are the unit normals in the x, y and 2 (depth) directipns, 
respectively. 



It should be noted that the solution -being presented assu-nes that the illuninated 
surface is in a sniall section (ySA << rj_) of the Fraunhofer region of the source 
field and that the field position is in the Fraunhofer region of the illuminated 
area (1*2 » r2 >> AA/x ), This point is emphasized by our A notation for the 
illuninated area. 

The mean square scattered pressure, which in general includes both 
coherent and incoherent sound intensity, is formulated: 



■=2*>= SSSS [V . K^v 



3 



X yx'y' 



(dxdx'dydy*) 
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where 



K 






u 


= X* - X 


[ 

r 


V 


= y' “ y 


i 

i 


^Pl 


= A?2 Ap^* 

plane plane 


{ 


* 


= complex conjugate 





( ) = ensemble average 

The ensemble average, < > , can be easily evaluated for our surface 
with its two-dimensional Gaussian distribution: 

, iK„ (C- C')v - Cka (cos 0 + cos 6-)] ^ Cl - C] 

\eZ /=e 1 ^ 

Define the ACOUSTICAL ROUGHNESS OF THE SEA SURFACE = g 

g = [kc (cos 0j^ + cos 6^)2 



O 
"5 ^ 



( 12 ) 



so that the mean scattered intensity can be written 



<Ap Ap *> = I 

2 2 pi (^)2 






X y X y 



i(K u + K v) , , 

e x y dxdx dydy 



If oiu* interests viere to be restricted to low roughness (g < 1) at this 
point, it would be useful to use the approximation e — (e“^) (l+gC), 

This would permit immediate separation into a coherent term, and an 

incoherent term which the integration would transform into tlie two dimensional 
surface spectral description, eqn« (^t) j (Weiner-Khinchine t'ransf ormation) . 



9 



The rarity of experimental measurements of two-dimensional surface spectral 
descriptions turns us away from this approach. 

0/ Instead, for scatter in the specular direction, it is convenient to 
try a series approximation :to the assumed isotropic Gaussian spatial correlation 

5,6 

function (Eq. 8). The integration then produces the relative scattered 
Intensity, (compared with the plane surface reflected intensity) specialized 
to specular scatter, in which 0 =0 a ^ q s F = 1 , g - 3^2 

123 ^ 

where = =i grazing angle. 



<pp*> 



<Ap2Ap2*> 

^plane 





(13) 



The two terns of Eq. ( 13 ) represent coherent scattered sound (the 
term) and incoherent scattering (the second term). A study of the tern 
in parentheses, defined as 

• S(g) 

shows that S (g) g 

and S (g ) 

L 

Since it can be assumed that L < AA, we can conclude that, in the limits, 
the specular scatter for a Gaussian sea will be: 

"• p* 

For Small Acoustical Roughness, g << 1 : <pp*) = e ' (l 5 ) 

2 

For Large Acoustical Roughness , g » 1 : (pp*) = ~^ = — (l6) 

;:k E Ai sin^t. 

In practice the coherent tom is dominant approximately up to g— 1 , and 

sometimes beyond that roughness, depending or the beam width, .Fig. 2 
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•mCf 00 

e ^ 2 

m=l 



m 



mmj 



for cr « 1 



for g » 1 



(lii) 



shows the relative behavlori- as a function of acoustical roughness. 

It is seen that for small acoustical roughness, the intensity of specular 

scatter is the factor times the miirror reflection , and the concept 

of a reflection coefficient is useful. 

However, it is important to observe also that the value of dB loss per 

bovmce is in general an inadequate, and in fact misleading, concept for the 

rough surface case. When g>^l the specific value of the scattering area 

enters into the asynptotic value of (pp*) which will therefore change from 

one bounce to the next as the sound wave continues to diverge and scatter, 

V/hen the Fraunhofer conditions are violated (large beam width) it is 

necessary either a) to add the scatter contributions from the several areas 

A A or b) to re-sol ve the problem, assiming a diverging wave at the 

T 8 

very beginning of the derivation, * The latter attack and its verification 

O 

by laboratory model° show that for specular scatter from an 

exponentially correlated surface there is virtually rx> difference between the 

Fraunhofer and Fresnel solutions if g<l. However, for g>l the scattered 

intensity may be significantly greater than predicted by the single Fraunhofer 

solution, T.iis is because, commonly, significant signal will exist over an area 

containing more than one A A as defined in our derivation. 

Laboratory models have been used in order to supplement the capabilities 
9 10 

of current theory, * We have recently studied the behavior of the coherent 
component for cross-wind, up and down wind specular scatter under widely-varying 
laboratory seas. Our Tseas" are produced by combinations of a paddle and five 
fans blowing air over the sxirface of our 25* x 6* x 8' deep anechoic tank. 

In this way we obtain water waves that give us scaled replicas of cesnpound sea 
and swell surfaces that are realistically time-varying, anisotropic, and near- 
Gaussian. The probability density function, spatial and tem-^oral correlations 



■**The Gaussian correlated surface would give similar results. 



and frequency spectrum compare closely to ocean measurements , except for 

scale factors for frequency, time and space lag, of course. 

Figure 3 shows three of the laboratory sea and swell spectra. The 

probability density functions of the wave height for the three cases are shown 

in Fig. U; the near-Gaussian probability density function is very close to 

that observed at sea.^^ Equation ( 13 ) states that the ratio of the intensity 

of the coherent component to the intensity reflected from a mirror-like 

surface should be e“^. Figure 5 sho;^s this ratio divided by e”S over four 

decades of acoustical roughness j for g^ 1 the ordinate is essentially unity. 

This excellent verification of theory, >^en tested \inder widely varying 

laboratory seas, has encouraged us to apply it to the real ocean. Other 

11 12 

laboratory model results ^ suggest that for g<l, at near grazing incidence, 
■a shadowing correction is not necessary. 



IV. EFFECTS OF NEAR-SURFACE BUBBLES 



Bubbles near the sea surface can modify the forward scatter of sound 
In two ways: (a) By near-surface refraction; the presence of all bubbles 

(resonant or not) will create a below-surf ace medium through which the speed 
of sound differs from that in the non-bubbly water; this change of speed will 
refract incoming rays, modify the angle of incidence at the rough-scattering 
surface and change the distance betweeen bounces in channel propagation. 

(b) By resonant absorption and scatter; resonant near-surface bubbles 
can cause attenuation that is generally unresolvable from rough surface 
scatter loss. 

The objectives of the present work are arbitrarily restricted to 
"surface*^ losses at frequencies below 16 kHz. Nevertheless, because bubble 
refraction depends on the total volume of air in bubble form, it is necessary 
to consider the numbers of near-surface bubbles of ail radii. 

In suite of the ^eat interest in the subject for aLmost 30 years, 
the literature of ocean bubble data is very thin: In-situ measurements of 

bubbles in off-shore breaking waves have been reported^^; bubble generation 
and persistance by an artifical wind over a laboratory tarJ< filled with sea 
water has been described^; inferred bubble densities have been calculated 
from certain AMOS data j in-situ acoustical measurements in 60 ft of iso- 
thermal water in Mission Bay, San Diego have been summarized^^j a laboratory 

study of surface bubbles has yielded data produced by colliding waves of 
17 

sea water© It will be necessary for us to extrapolate beyond these bare 
measurements available. 

A. Dependence on Radius 

Because of experimental difficulties, in no case are there published 
jji-situ measurements of ocean bubbles of radius .greater than 200 
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microns (resonant frequency less than approximately 15 kHz) or less than 
roughly 20 microns (resonant frequency greater than 150 kHz). Most of the 
available data fit approximately the solid line in Fig. 6 > plotted on 
the basis of fractional volume of air in bubble form to water u(R)dR. From 
the fractional volume it is easy to calculate the number of bubbles per unit 
volume, in radius increment dR, by using 



n(R) dR = 



u(R) dR 

4/3nR^ 



( 17 ) 



There is some evidence that the bubble populations break up roughly into two 
groups: those smaller than approximately 60 micron radius postulated as 

attached to motes and due principally to biological activity, cosmic rays, 
captured aerosols and photosynthesis of phytoplankton, and the larger bubbles 
which we assume come from breaking waves (and decaying matter on the sea 
floor in shallow water) . 

,v- To calculate the near-surface refraction due to bubbles we have had to 
extrapolate (dashed lines) the data~based curves of Fig, 6, which were 

measured with incipient breaking waves at wind speeds U-6 knots at 3,3 meter 
depth. For radii smaller than 30n we extrapolate to one micron radius at the 



KoR 



- 1/2 



dependence that has been found for the smallest measured ocean 



bubbles^^ and in tap water^^j below 1' micron, at slope following behavior 

19 20 

loosely suggested by laboratory cavitation studies, * For bubbles larger than 
60u we have used the dependence observed^ to 100 microns and have 



lU 
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_3 

extrapolated to largest radii at Blanchard and Woodcock's slope of K^R 

The radius at which one shifts to the suggested slope for largest bubbles 

is, admittedly, somewhat arbitrary. That slope starts at about 50 microns 

in the breaking wave data; however, a re-analysis of Glotov's data shows that 
-3 

the K^R behavior is confirmed for R > 120p . The constants for the slopes 
of Fig, 6 are shown in Table 1. 

TABLE ONE 

Estimated Fractional Bubble Air to Water at 3.3M depth for Winds Approximately 
\ 6 Knots 



Radius Range 
Microns 



Dependence 

u(Rq,3.3M) 



Constant [dimensions] / u(Rq) dR 











0 








< 


1.0 


K,R« 


K. - 




0 






1 0 


1 


1.0 <_ 


R 


< 


30. 


K-R 


K., = 




0 






2 0 


2 


30. _< 


R 


< 


60. 




Ko = 




0 






3 0 


3 


60. 




< 


100. 




K/ = 




0 






4 0 




100. <_ 








K R 


K^ = 




0 






5 0 


5 



= 1 V TO 






-9 

1 t; 10 



1.5 X 10"^ 13.8 X lO"^ 



8.5 X lO"^ [y°] 5.8 x lO”^ 



-12 -2 -9 

2.5 X 10 [\i ] 7.7 X 10 



-4 +2 -9 

2.5 X 10 [y ] 12.4 x 10 



= 4.1 X 10 
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Total Effective U 



B, Dependence on Depth 

For our purposes, interest lies in the bubble densities (and fractions) 

as a function of depth at a fixed frequency, rather than for fixed radius. 

Studies in well-mixed ocean water suggest that for bubbles larger 

than approximately 30 microns (sea level resonant frequency less than 100 kHz) 

buoyant forces and turbulent entrainment determine the bubble life time 

rather than gas diffusion through the bubble walls. If we formally make the 

simplifying assumption that these bubbles can be lost only at the ocean surface, 

then bubbles of surface (3.3 Meter depth) resonant frequency f^, radius R^, 

will be found at depth Z in the same numbers, but at the isothermally* 

1/3 

compressed radius, resonant frequency f^^ where 

P “ 1.33/(1 + O.IZ). 



For resonant air bubbles 1 >;Rq=— 
y = 1.4, P = 10^ ^1+2/10 dyne/cm^, °f> « 
meters , 



where 1.5x10^ cm/sec 
1 . 0 gm/ cm , Z = depth in 



The bubble fraction at depth Z will be 
u(Roz.Z) = n(RQ2,Z)(4/3 



Analytically, our assumption for large bubbles is: from the surface 

to depths where the entrained bubble drag force is equal to or greater than 
the buoyant force 

n(R^,3.3M) dR = nCR^^.Z) dR ^ 



*The 

restriction to 
changes will cause 



isothermal" waters is not severe, 
negligible volume changes. 



since small temperature 
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80 that following isothermal compression, 

“ [u(Rq,3.3M)][6] . (^) 

These bubbles, however, have a resonant frequency which is not f^, but, 
because of changed radius and pressure, is 



OZ 



^0 ^ 



-5/6 



' We can now calculate the gas fraction of bubbles at the desired 
resonant frequency f^ at depth Z, u(f^,Z) , by assuming that the frequency 
dependence within a given band of bubble sizes is the same as at the surface, 
that is. 



u(f^,3.3M) 



k * f.” and u(f._,Z) 
n‘>-^ 0 * 




f 



n 



(JZ 



where n is the frequency dependence within a particular frequency band. 
Therefore 



u(fQ,Z) = 



n 






(19) 



Finally, from (18) and (19) 



u(fQ,Z) = 



(^5/6 n+1^ 



u(fQ,3.3M) 



and, taking account of the change in radius. 



n(f^,Z) 



(g5/6 n + 5/2^ n(fQ,3.3M) 



( 20 ) 



( 21 ) 
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The predicted depth dependences for the three most important bands 
of surface radii are shown in Table 2# 

TABLE 2 

Fractional Air in Bubble Form and Number of Bubbles per Unit Volume 
as a Function of Depth for Large Bubbles 



^0 

t 

microns 


^0 

kHz 


n 


u(fQ,Z)/u(fQ,3.3M) 
p5/6 n+1 

m 


N(fQ,Z)/N(fQ,3.3M) 
-5/6 n + 5/2 

“ P 


30 to 60 


60 to 100 


+1 


jll/6 


jlO/3 


60 to 100 


32 to 60 ■ 


-1 






> 100 


less than 32 


+3 







It is clear that the assumptions that lead to the values in Table 2 
can only be considered as a perturbation analysis of the change of bubble 
populations with depth since they imply inconsistencies, particularly at the 
edges of the bands, as depth is increased. Nevertheless, the only bubble 



density data available 



plotted against the absolute pressure, 1.33B 

+2 



for four frequencies between 40 and 96 kHz show slopes of approximately B , 

1.7 3.3 

in reasonable agreement with the predictions of B * to B * . 

When we plot the frequency band contributions to the total bubble fraction, 
( / u R dr), as a function of total pressure parameter B the total bubble 

fraction U(Z) at wind speed 6 knots is found to vary (with less than deviation) 

+2 +1 

as B from the surface to approximately lOM depth, blending into B 

at greater depths (near 20 Meters) . Therefore the dpeth dependence that we 

will use for the near-surface refraction calculation is: 



= 4.1 X 10“® 
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U(z) 



( 22 ) 



C. Dependence on Wind Speed 

It has been found ^ that the number of bubbles shows "an abrupt 

Increase ... approximately ... exponential" for winds above about 6 knots. 

21 

The comment agrees roughly with the curve of foam coverage vs. wind speed 
which can be written approximately as: 

U(Z=0,W) = U(Z=0,6 knots) [1 + 10"^(W/10)^ (23) 



where W = wind speed in knots. 

There is a second wind effect: the larger, wind-generated bubbles will be 

brought to depths by increased turbulent entrainments . A complete knowledge of 

this behavior may be possible when there is a better understanding of under- 

22 

water motion as a function of wind speed, insolation and depth. 

However, the factor of bubble entrainment by water motion can be crudely 
evaluated by the following simple argument: 

A bubble will be drawn down from the surface when the drag force is 
greater than the buoyant force. 

The co.XLTective velocity in the medium will probably be Gaussian distributed. 
Nevertheless, the mean downward convective velocity will increase with wind 
speed as does the velocity spectrum and the depth of the well-mixed layer. 

We assume that the mean convective velocity is proportional to the wind 
speed and that the mean number of bubbles drawn downward under a given 
condition of underwater turbulence is proportional to the mean velocity; the 
dependence on radius will be ignored at this time. This "zero order" 
estimation of the effect of turbulent entrainment, combined with eq.(23) 
for the increased foam available leads to the wind dependence: 



U(Z,R ,W) 

=(W/6) [1 

U(Z,Rq, 6 knots) 



+ 10"^W/10)^] 



( 21 ;) 
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D. Summary of Bubble Dependences 

Combining the dependences on radius of bubble, bubble depth and wind 
speed, we get the total bubble fraction: 

U(Z,W) - [U.l X 10"®^^] [W/6J [1 + lO'^(WAO)^] (25) 

It is appropriate to point out that other bubble factors exist in 
addition to the effects that we have considered. Two significant elements 
are the dependence on dayli^t and precipitation; others are dependence on 
water temperature and gradient, percentage gas sat\iration, marine biological 
activity. As examples: Blanchard and Woodcock have shown that snow and 

rain produce very large numbers of bubbles. Their penetration in laboratory 
experiments was no more than 2 cm below the water surface but it would be 
expected that these bubbles would penetrate much farther in the presence of 
convection cells accompanying storms at sea. Small bubbles (less than about 
60 V ) are about twice as populous in the daytin.e as at night. It has been 
postulated that the increased numbers may be due to photosynthetic action, 
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E. Theory of Dispersion and Absorption by Bubbles 

Following classical developments an incident acoustic pressure, 
of wave length very large compared to the bubble radius R, having a velocity 
potential ij)^e will produce a total field , 



N 



<f>(ir) = 0. (ir) + X (~r) 



ik 1 ir “ ir I 
e n 



( 26 ) 



where ir = coordinate position of scattering region 



or-ir = 



n 



(0 






th* 

distance between' n bubble and center of scattering region 
2n/X = propagation constant in bubble-free water 

In 

(woA)^ - 1 - i6 

bubble radius 

incident sound frequency 

resonant frequency of bubble 

•* • 

damping constant of bubble 



The average over all possible sizes and configurations of bubbles is 

<m> « Or) *///g (g 



n 



where 



(27) 



G (ir ) = y* Cn (ir,R)] [" 1 

R LC‘^o(R)A 3 -l-i6(u;,R)J 
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(U) 



Equation (27)is the integral equation formulation of the differential 



equation 



+k^( 0 ) = “ 4 irG< 0 ) or 



<^(|> + fc ^ (r) <(||> = 0 



( 28 ) 



where 



kj^(r) = k 



ri4-^>-k+f • f 

L k^ J . ^ -^R KA') - 1 -i 



6(R) 



The quantity, — : = 1 + — — , is the index of refraction in the bubbly 

iC ^ ^ 

region. Because of the presence of the damping term, the index will have 
an imaginary part which comprises the attenuation due to absorption and 
scatter in the bubbly water in addition to the real part, defining the local 
speed of sound. These terms are: 



Re; {k. (r)| = k + 



27T 






R n(R)[(i^)^-l]dR 



(29) 



Im: (r)} = ^ / 

R 



6 R n (R) d R 



[e) -^] 



(30) 



+ 6 
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F. Predictions of Refraction by Bubbles 

I 

The dispersion relation Eq. (2?) > since it predicts a decrease ^ 

of sound speed caused by bubbles of resonant frequency greater than the 
Incident frequency and an increase of speed for bubbles with lower resonant 
frequencies, would require precise information of bubble populations for 
an evaluation for all sound frequencies- However, if our Interest is restricted 
to sound frequencies less than 16 kHz, integration of Fig. 6 shows that 
bubbles under that resonant radius (R^ < 200. microns) represent 

approximately 93% of the total fractional air volume. For the calculation j 

of refraction at frequencies less than 16 kHz we consider only the effect of 

bubbles of resonant frequency greater than the driving frequency# [ 

' 2 

When >> Lj , we drop 6 in the denominator of Eq. (29i,use 
definition of u (R) to obtain the simplified form: 






iO << W, 



(31) 



Since the expression in parentheses is nearly unity, we find 



c 1 - 8.0 X 10^ U (l + Z/10)“ 



‘] 



(32) 



where U(Z,U) is given by equation (25) • 



i 

I 
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The significance of this departure of the speed from the bubble-free 
value is evident when we insert Equation (25) and look at the gradient 

of the near-surface propagation speed. For the reference wind speed of 

-4 

6 knots the gradient of the velocity, proportioanl to (1 + O.IZ) ^ has a 
surface value of 0.26 sec decreasing to 0.016 sec ^ by 10 Meter depth. 

Since the velocity gradient due to pressure in isothermal water is 0.018 sec**^ 
our conclusion must be that ray curvature in approximately the top 10 meters 
is very strongly influenced by bubbles even for winds as low as 6 knots. 

Because of bubbles, rays approaching the surface will refract upward, 
meeting the surface in a cusplike ray path, and making more bounces for a 
given horizontal range than for a bubble-free region. 



G« Predictions of Attenuation due to Bubbles 

The sound pressure attenuation formula for the bubbly water can be 
written as 



^'-Ofx i(cjt“k, x) 
*- e e b 



where 






Using Eq.(30)and rewriting in decibel notation, we get the attenuation 



a, = 8.68 a 

b 



(8.68) W J 



“ (6/6^) R“ n(R,Z,W) dR 






+ 6 



(33) 



where 6^ = radiation damping constant = k 



ta 



1.36 X lO"^ ^ 1+Z/lO for air bubbles 
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Integration yields 



= [8.687tV\(R ,Z,W)]/6 

D 0 O 1 



which can be re-expressed as 



>- 1/2 



a, = 1.50 X 10^ (1+Z/lO) u(f.,Z,W) db/Meter (3U) 



The conclusion of the calculation of "surface" loss in db/bounce 
due to bubbles requirer integration of Eq. (3U ) over the surface duct ray 
path of varying curvature and varying attenuation: 






where L is the depth at which the ray direction is horizontal. Inserting 
a^ and 



ds = 



dZ 



sin 0 



r>J 



ftdZ 



. 1 



2 j2ft(L-Z) ^ V 2(L-Z)G 



dZ 



where ^ = local radius of curvature, C-(Z) is the total local gradient, and 
c = local velocity, M/sec, we obtain 






= 1.06 X 10' 



0 u(f^,Z,W) 



dZ 



^ (L-Z)^'^" (1 + O.IZ)^'^^ G^'^^(Z) ‘ 



( 35 ) 



For frequencies below 16 kHz, the expression for the bubble fraction, 
from Table 2 and Section C, is 



u(fQ,Z,W) = 0.8 X 10’^^ f^^ + 10"^(W/10)2j 



I where W is in knots and fQ is in Hz, 
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V. SURFACE DUCT PROPAGATION - COMPARISON WITH EXPERIMEin* 

1 

Bucker has recently solved the surface duct propagation problem by 
a wave theory approach. His modal components, which correspond to up and 
down propagating waves, with appropriate ray directions, lose energy when 
they interact with the surface and by leakage through the thermocline. Our 
proposal here is to introduce the results of our development in Sections III 
and IV into Bucker 's solution. Two changes are necessary: the velocity 

gradient in the duct will now be due to bubbles as well as temperature and 
pressxire effects; the loss at the surface will l>e made up of the losses per 
bounce due to below-surface extinction of sound by bubbles and simple, at-surface, 
coherent specular scatter (because, in general, g<l), calculated for each 
contributing mode. The surface "reflection coefficient" is thereby directly 
attributable to the wind velocity which is the sole source of the rms wave 
height* in the examples chosen and a principal cause of the bubbles. 

First, let us look at how the losses vary with wind speed and ray angle 
of incidence. Our loss calculations are based on eqs. 1^, 35> 36. Our ray s 

r 

paths are calculated from the bubble-caused gradient using eqs. 25, 32 together 

with a temperature plus pressure-caused gradient of O.OliS sec”^. Figure 7 

illustrates the relative contributions of surface loss and integrated bubble 

loss for a range of winds with a ray that is horizontal at UO M, ^ = U. 38 *. 

below I 

The frequency is'^ 16 kHz. It is observed that the bubble contribution | 

I 

dominates only for low wind speeds. \ 

j 

*In the data report^"^ the "seaman's eye" estimates of the peak-to-peak wave * 

height, H, are in good (805o) arrcement with the calculation from • equation (6) | 

when it is assumed that H = 6u- , for Gaussian wave systems. However, our ; 

laboratory experiments suggest that, in general, the total height due to s 

swell and "sea" should be used to determine o*. j 

I 

! 
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Figure 8 shows that the bubble loss per bounce is not very sensitive 
to the depth of the horizontal ray (or angle of incidence at the surface); 
the loss per bounce at the surface, however, increases as the square of the 
jine of the grazing angle. 

' For higher frequencies it would be expected that the neglect of the 
Incoherent surface scattering term would no longer be possible; the forvrard 
scatter loss would reach an asymptotic value which would be smaller than the increased 
bubble loss. For frequencies lower than o\ir example the coherent scatter 
approximation, e“®, becomes even better, and would be applicable to higher 
sea states; the bubble loss would be expected to continue to be smaller than 
the surface loss. 

The experimental propagation data for our study are again from Runs 73 
and 99 of the operational tests of the SQS-26 sonar^"^ in which both source 
and receiver were in the duct. The wave heights are such that g<l for ail 
affective modes of the duct solution. 

Since Bucker's duct solution assumes constant velocity gradients, we 
lave had to approximate the calculated total gradient, which goes as^“^ , 
ly splitting it into constant sections and G 2 (Fig. 9). The near-surface 
section, G]_ , has been given a value twice the uncorrected gradient due to 
temperature and pressure and the remainder, G 2 , has been assumed to have no 
subble influence. The values are given in Table 3* In this way, the duct 
sropagation loss cui*ves have been recalculated by Bucker. They compare well 
rith the data (Figs. 10 , 11 ). The dashed lines in Figures 10 and 11 are the 
ituch poorer predictions based on the equations of project 
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TABLE 3 



Bcfba Tor Two Surface Duct Experiments 



Run No. 




Vs 


H °i 


®2 


W 






M/sec 


M sec"^ 


sec“l 


knots 


73 




1513.1 


36.6 0.096 


0 . 0 U 8 


17 


99 




1510.9 


51.8 0.102 


O.O 3 I 1 


15 


where 




= velocity 


at surface 








Zt 


* depth of 


therraocline 










* velocity 


gradient, surface to 


8 yards 








" velocity 


gradient, 8 yards to 


thermocline 





Other, higher values of the assumed bilinear gradient and a closer- 
fitting tri-linear gradient have also been considered. However, the resulting 
propagation prediction curves show excessive losses compared to experiment# 
Since, at these wind speeds, the bubble loss is less than the surface loss, 
the velocity gradient can be a major influence. This suggests that our total 
bubble populations, U(fQ,Z,W), may be too high, perhaps by a factor of two, 
possibly due to the wrong wind speed dependence having been assumed. 
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VI. CONCLUSION 



The purpose of this paper has been to show that for a large range of 
ocean surface roughness, g<lj it is possible to predict surface duct 
propagation by using the simplest surface scattering theory (coherent term 
only) and estimates of bubble presence as inputs to the Bucker modal solution. 

■ Nevertheless, it is the necessity of kncv?ing more about the ocean 
parameters of sound propagation that has been demonstrated by this paper. 

It is hoped that the present work will serve as a spur for more fundamental 
studies to replace the speculative assumptions that we have been forced to 
make about bubble dependence on wind speed. 
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